Diffusive transport of light in two-dimensional granular materials 
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We study photon diffusion in a two-dimensional random packing of monodisperse disks as a simple 
model of granular material. We apply ray optics approximation to set up a persistent random walk 
for the photons. We employ Fresnel's intensity reflectance with its rich dependence on the incidence 
angle and polarization state of the light. We present an analytic expression for the transport-mean- 
free path I* in terms of the refractive indices of grains and host medium, grain radius, and packing 
fraction. We perform numerical simulations to examine our analytical result. 
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I. INTRODUCTION 



Pebbles on a sea shore, sand, rice, and sugar, are a few 
examples of ubiquitous granular systems [I|, [2[ . Granu- 
lar media consists of discrete particles of size larger than 
100 /ini, interacting with each other through dissipative 
contact forces. In the absence of an external drive, par- 
ticles rapidly lose their kinetic energy. However, granu- 
lar materials under external forces exhibit transition be- 
tween a liquid-like and a solid-like state. In order to 
understand jamming transition |3j, size segregation [4J, 
convection rolls, pattern formation and dynamical insta- 
bilities [5] of granular systems, probing the micron-scale 
dynamics of constituent particles is important. Although 
granular systems are opaque, diffusing wave spectro scop y 
(DWS) [6] non-invasively probes their dynamics |7H22|. 

In a turbid material, light experiences many scattering 
events before leaving the sample, and the transport of 
light energy is diffusive [23|. Consequently, the photon 
can be considered as a random walker. The structural 
details of the opaque medium reflect in the transport- 
mean- free path /*, over which the photon direction be- 
comes randomized. Moreover, the dynamics of scatterers 
leads to the temporal intensity fluctuations in the speckle 
field of the multiply scattered light. Utilizing the inten- 
sity auto-correlation, DWS determines I* and the mean- 
squared displacement of the scatterers. DWS has been 
used to study colloidal dispersions B , liquid crystals [24j , 
biopolymers 25j , and foams [26|, |27| ■ 

Multiple light scattering from grains has been in- 
voked to reveal their relative motion. Three-dimensional 
gravity-driven granular flow [7( , gas-fluidized beds [8j , 
water- fluidized beds jlljl , vibro-fluidized systems [9, Il6r - 
|21| . avalanche flow [l(|, creeping motion |12J, dilation 
due to temperature variations [14] . and response to a 
localized compression force [151 ] have been thoroughly 
studied. As a key parameter, I* was measured. For 
glass beads dispersed in water, I* « \AR — 1QR for 
80 /zm ^ R ^ 200 /an. These samples had a pack- 



*Electronic address: miri@iasbs.ac.ir 



ing fraction <fi w 0.64 |22j. For glass spheres of radius 
R = 47.5 urn dispersed in air, I* « 15i? is reported Q- 
Crassous [13| developed a ray-tracing program to access 
I*, but only for packing fraction 4> ~ 0.64. 

It is natural to study I* as a function of the refractive 
index of grains rij n , refractive index of the host medium 
n ou t, average grain radius R, and packing fraction <fi. In 
one approach, expansion of electromagnetic fields in a 
series of vector multipole fields [28| |. or other accurate 
techniques can be invoked to simulate the speckle pat- 
tern. To extract I*, simulations must be repeated for a 
large number of realizations of the random system. Nu- 
merical calculations are quite demanding as each sample 
contains hundreds of grains. In another approach, we 
focus on the elucidation of mechanisms underlying the 
random walk of photons. As grains are much larger than 
the wavelength of light, we rely on ray optics approxima- 
tion. Within this framework, the role of reflection and 
total internal reflection phenomena in abrupt change of 
photons' paths, can be pictured. 

We study photon diffusion in a two-dimensional ran- 
dom packing of monodisperse disks as a simple model of 
granular material. We assume that the grains are ho- 
mogenous and transparent. We employ ray optics ap- 
proximation to follow a light beam or photon as it is 
reflected by the grains with a probability called the in- 
tensity reflectance. The photon's random walk based on 
the above rules is a persistent random walk 29]. This 
shows that the diffusive transport of light in granular 
media and foams [30, [31| are much similar. As an exten- 
sion of our previous study [31j , here we take into account 
that the intensity reflectance depends on the incidence 
angle and the polarization state of the light. Writing 
master equations to describe the photon transport, we 
obtain analytic expression for I* as a function of model 
parameters n,„, n outl R, and <p. We perform numerical 
simulations to examine our analytical result. 

Our article is organized as follows. In Section Ql] we 
introduce the model. Photon transport in a random 
packing of disks using Fresnel's intensity reflectance is 
discussed in Sec. IIIII Discussions, conclusions, and an 
outlook are presented in Sec. [IVJ 



II. MODEL 

Following a step-by-step approach to reality, we delib- 
erately focus on two-dimensional granular systems. Our 
model granular medium is a random packing of circular 
disks. All nonoverlapping disks have the same radius R, 
and cover a fraction <f> of the plane. The refractive indices 
of grains and host medium are rn n and n ou t , respectively. 
We assume that n in > n out . 

Grains (> 100 /im) are larger than the visible light 
wavelength (400 — 700 nm). Therefore, we employ ray 
optics. A light beam or photon experiences transmission 
or reflection as it hits the surface of a grain. We denote by 
r _j.i the intensity reflectance for photons moving in the 
host medium and hitting a grain. Similarly, we denote 
by Ti^ the intensity reflectance for photons moving in 
a gain and hitting its surface. According to the Fresnel's 
formulas, both r ^i and r,_,. depend on the polarization 
state of the light, incidence angle 7, n out and rii n , see 
Appendix [X] Indeed the Fresnel's intensity reflectance 
n^o(7) is 1 for 7 > 7c, where 



+ 7T + 27 



7 C = arcsin(n „ t /ni„) 



(1) 



is the critical angle. Our model does respect the total 
internal reflection phenomena. 

It is instructive to consider a toy model, a one- 
dimensional lattice of (point) grains. On hitting a grain, 
a photon will be either reflected by probability r or per- 
sist on its direction of motion [32J. This leads to a per- 
sistent random walk of the photons, where the walker 
remembers its direction from the previous step [29J, |33| . 
Now it is clear that implementing the rules of ray optics 
for photon diffusion in two- and three-dimensional granu- 
lar media results in a generalized persistent random walk. 
Firstly introduced by Furth as a model for diffusion in 
a number of biological problems 34] , and shortly after 
by Taylor in the analysis of turbulent diffusion [35| , the 
persistent random walks are employed in polymers |36l |. 
Landauer diffusion coefficient for a one-dimensional solid 
37| , and general transport mechanisms [38( . 



III. PHOTON TRANSPORT IN A 
TWO-DIMENSIONAL GRANULAR MATERIAL 

A. Analytical treatment 

To simplify our analytical treatment of photons ran- 
dom walk, we further assume that (i) The length of pho- 
ton steps outside (inside) the grains is L out (Lin), i-e. 
we neglect the fact that the length of photon steps are 
not equal, (ii) A photon which transmits into or out of a 
grain, does not change its direction of motion. In other 
words, we neglect the fact that the angle of refraction is 
not equal to the angle of incidence. 

In the two-dimensional space, each photon step can be 
specified by an angle relative to the x axis. Consider a 




FIG. 1: (a) Path of a photon moving in the host medium 
and hitting a grain with an incidence angle 7. (b) Path of 
a photon moving in a grain and hitting its surface with an 
incidence angle 7. The step length inside the grain is 2R cos 7, 
where R is radius of the grain. 



photon which moves in the host medium along the di- 
rection 9 + -k + 27 and hits a grain with an incidence 
angle 7, see Fig. [TJa). The photon will be either re- 
flected to the direction 9 by probability r ^i (7) , or enter 
the grain. The probability distribution of the random 
variable < 7 < 7r/2 is [31j 



Foutil) = COS7. 



(2) 



Similarly, a photon moving in a grain along the direction 
8+tt+2"/ and hitting its boundary with an angle 7, will be 
either reflected to the direction 9 by probability r^ (7) , 
or leave the grain along the initial direction 6 + n + 27. 
Now the probability distribution of the incidence angle is 



*i»(7) 



cos-f/sin-fc |7| < 7 C 







otherwise, 



(3) 



see Appendix iBl Note that rii n > n out and ri_> (7) = 1 
for 7 > 7c. Thus 7 < 7 C ensures that a photon moving 
in the grain is not trapped forever. 

We let P° ut (x\9)dxdy (P™ (x\9)dxdy) denote the prob- 
ability that a photon moving along the direction 9, ar- 
rives at vicinity dxdy of position x = (x, y) after its nth 
step outside (inside) a grain. We express the evolution of 
P° ut (x\9) and P n n (x\9) by following master equations: 

*S+i(x|0) = \ J I P° n ut {*-L out ee\9 + it + 2j)F mt (7) 
xr o ^,-(7)d7 + ti-^ o i^ n (x-£ <m tee|0), (4) 
P£+d*\0) = \\\ ^"(x - L m e e \9 + t: + 2 7 )P m ( 7 ) 

xr<_> (7)d7 +i o -*i*T l *(x-L in efl|0), (5) 

where eg — (cos 9, sin 9) is the unit vector along the direc- 
tion 8, L out (L in ) denotes the average length of photon 
steps outside (inside) the grains, and 



(1 -r -n{')))F mi .t{i)dri, 



ti_>o = \ J C (1 - r^ ( 7 )) F ln { 1 )d 1 . (6) 



The first term on the right-hand side of Eq. §4$ represents 
the reflection of the photon with a probability r _>j(7). 
On arriving at position (x—L out cos 9, y—L out sin 9) along 
the direction 9 + ir + 27, the photon changes its direction 
by an angle ir + 2j according to the probability distribu- 
tion F ou t (7) . The second term focuses on photons which 
move inside the grain and transmit to the host medium 
with a probability ij_j. - The photon which has arrived at 
position (x — L out cos 9, y — L out sin 9) along the direction 
9, moves ballistically to position (x,y). Equation ([5]) has 
a similar interpretation. 

We are interested in the second moment of photon dis- 
tribution with respect to the spatial coordinates x and 
y. To evaluate the moments of an arbitrary distribution 
function P n (x. y\9), we utilize its associated characteris- 
tic function 12911 



P n (u x ,co y \m) 
Indeed 

(x kl y k2 ) n = 



:x -P n (x,y\9)dxdyd9. 



(7) 



H) 



x kl y k2 P n (x,y\9)dxdyd9 
fci+fc3 afci+fc2 Pn( j| m = ) 



du&ddfr 



,(8) 



where k\ and ki are either zero or positive integers, and 
u = (u x ,u y ). 

The Fourier transform of master equations (|4]) and §5§ 
are 

P°fi(«,a|m) = 

+ OO 

£ i k e-' ka J k (uL out ) [ii_, P™(u;, a\k + m) 

k— — oo 

+ c£? fc P£ rt (w J a|*+m)' 
P^ 1 («,a|m) = 

J2 i k e- tka J k {ioL m ) [i ^P° ut (u;, a\k + m) 



k=—oo 



•c£ ifc P*>,a|fc+m) , 



(9) 



where uj and a are the polar representation of the vector 
w = (cu x ,ujy), Jk is the fcth-order Bessel function, and 



(-1) 



m-\- h 



~"m,k 



J-i) m+k f c _ i{ 



e-^ m+2k ^F 0Ut { 1 )r ^ l { 1 )d ll 
e- l{2m+2k ^F m { 1 )r^ { 1 )d 1 .{\Q) 



We introduce the Taylor expansions 



PSi(«.«H 



ra+1 



(u, a\m) 



Q°ot(a\m)+iu;L out Q™t(a\m) 

, , 2 T 2 

2 t *2,n "r ••! 

Q™ n {a\m) + iuiL in Q™ n (a\m) 
_w^, Q i„ + _ (u) 



Now we insert Eq. (fTT|) into Eq. (|9]). Using the Taylor 
expansion of the Bessel functions and collecting all terms 
with the same power in w, we obtain a complicated set 
of recursion relations for Q° u *{a\m) and Q™ n {a\m), see 
Appendix [C] There is an elegant method to transform 
this set of coupled linear difference equations to a set 
of algebraic equations: the z transform [29|,[39|]. The z 
transform Q(z) of a function Q n of a discrete variable 
n = 0, 1, 2, ... is defined by 



q(z) = J2q, 



n=0 



(12) 



The z transform of Eqs. (|C1|) - (|C3[) leads to a set 
of algebraic equations whose solutions Qf l (z\a, m) and 
Q° ut (z\a, m) are reported in Appendix [Cl 

From Eqs. © and (fTTj) it follows that the first and 
second moments of photon distribution are 

(x) n = L out Qt%(0, 0) + L in Q^ n (0, 0), 



,7T 



, 71" 



(y)n = L out Ql%-^) + Lm Q\ n n (-,0) : 
(x 2 ) n = L 2 out Q°%(0, 0) + L^gj7 n (0, 0), 

,.2\ ._ r2 



ur>- = ^„/Q^(|,o)+iiQ™(|,o). (13) 



Powered by the analytical expressions for the inverse z 
transform of Qf l {z\a, m) and Q° ut (z\a, m), we find that 
in the limit n — > 00 

(»)n = (y)n = 0, 

(x 2 )n = {y 2 )n = {w out L 2 out + w ln L 2 m )n, (14) 
where 



Wout = 



to—¥i \H— >o 



to^ri t ^Z— »0 



+ 



^in-^out^l —xH—>o ~r h_ >o (W— ^o^o— M t Cq j \i ^0,1/J 



(1 — Cj^ )\~ c l.o) — to^riti^to 
i- J vn,'- J uV'0— : >%ki— : ro T C^ _j (ti— yoto— H T Cq I (1 — c o.l)J 



-0,0 

2 



V 1,0 s\ c l 0/ 0— HH— >o 



to— H i t'i^o 



(15) 



+ u 



-0,0 



in -L J out<'o—>iti—>o ~r C-y _-± \ti—> t ^i + ^0,1 [1 c o.l /J 



/I _ „out\(-\ _ „in \ _ f .y-. 
in Lout^X —\to—¥i "T" ^o— >i (^i— >o^o— >i H~ C 0,1V^ c l,0 /J 



V 



(l-c$)(l-c£ ) 



^o—^i^i- 



t, 



+ u 



(16) 



The task is now expressing the time nr spent for the 
n steps of the random walker. Indeed 



JoutTout i Jin^~ii 



(17) 



where f out (/j n = 1 — / ou t) is the fraction of time that 
the photons spend outside (inside) the grains, r out = 
n utL ut/c (Ti n = n in L m /c) is the average time spent 
to make a step outside (inside) the grains, and c is the 
velocity of light in vacuum. To evaluate /,„ , we consider 
a photon hitting a grain with an incidence angle 7. The 
probability of m internal steps before leaving the grain is 
ti^o{l)[i"i^o{l)] m -1 - Thus the photon spends a time 



m— 1 



mTinU^o (7) [rj-,.0 (7)]' 



i/u^oh) 



inside a grain before leaving it. Averaging with respect to 
the probability distribution Fi n (^), we find that a photon 
spends a time 



Tin-Tin \f ) 



- 7c u^o{i) 

inside the grain. It follows that 

Tout 



dj = T in (- ) 



(18) 



J out 
Jin 



7~in \ ^ . / \ 'Tout 

rJ(^) 

Tin\~ / ~r Tout 



(19) 



We also note that </> = Li n /(Li n + L out ). Figure [TJb) 
suggests that L in = (2i?cos7) where 7 is the incidence 
angle of photons moving in the grain, and here (...) de- 
notes averaging with respect to the probability distribu- 
tion F in (j). Hence we find 



Lout = i?(cos7 c + 7 C / sin 7 C )- 
L in = i?(cos7c + 7c/sin7c). 



(20) 



As already mentioned, in the long-time limit nr — > 00, 
the behavior of the mean-square displacements is purely 
diffusive, i.e., 



(X 2 )n = (y 2 )r 



2D™. 



(21) 



In two-dimensional systems, the transport-mean-free 
path is defined via 



I* = 2D/v mi 



(22) 



where v m is the transport velocity of light in the medium. 
To a first approximation 



(1 



n ou t 



(23) 



where c/n ou t (c/ni„) is the velocity of light in the host 
medium (grains) which cover a fraction 1 — ((f)) of the 
plane. Now we utilize Eqs. (TT4"l) - (j2"3")) to derive the diffu- 
sion constant and transport-mean-free path in the two- 
dimensional granular material 



D = 



- Rc\n in /n out aicsm(n out / 12^) + ^ l-(n ou t/ n in )'- 

/ 4-^s9 \ n out(^r) +n m ( T ] — } 
x (wout(—^) 2 + w m ) * ^f", 24 



R(n m /n ut&rcsm(n ou t/n in ) + y/ 1 - (n out / 'n in ) 2 
1 



nlut^ + nUl-h-) 



(!-</ 



(25) 



We emphasize that w ou t, Win and (jr^-) can be explicitly 
expressed in terms of physical parameters n ou t, nin, and 
<j>. The incident electric field can be either perpendicular 
_L or parallel || to the two-dimensional plane covered by 
the grains. The intensity reflectance depends on the po- 
larization state of light, thus the diffusion constants D± 
and Dm are not equal. 



B. Numerical simulations 

We carry out numerical simulations to examine our 
analytical result for the diffusion constant D. Using con- 
tact dynamics simulations [40| , we generate homogeneous 
and monodisperse random packing of disks. We study 
samples with packing fraction cf> £ [0.15, 0.25, ..., 0.65]. 
Each sample consists of 10 4 nonoverlapping disks. We 
let 10 4 photons perform a random walk in each sam- 
ple. We launch the photons in a direction specified 
by angle 60 ■ We repeat the simulation for all angles 
6»o £ [30°, 60°,..., 360°]. We implement Fresnel's formulas 
and Snell's law as a photon hits a grain surface. Follow- 
ing a standard Monte Carlo procedure, we generate the 
trajectory of each photon, and evaluate the statistics of 
the photon cloud at different times to access the diffusion 
constant D. For improving the speed of our ray tracing 
program, we adopt the cell index method commonly used 
in the molecular dynamics simulations, see Refs. [30| and 
[3l| for details. 

As an example, we consider the glass disks (n.j„ = 1.5) 
immersed in the air (n ou t = 1-0). Figures EJa) and[2jb) 
demonstrate D± and Du as a function of <fi, respectively. 
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FIG. 2: (Color online) The diffusion constant (a) D±, (b) D\\ 
(in units of the disk radius R times the velocity of light c) as 
a function of the packing fraction <f> for the case m n = 1.5 and 
riout — 1.0. (c) D± and (d) Du for the case rii n = 2.0 and 
riout — 1-34. 



As another example, we consider glass disks (n in = 2) 
immersed in the water (n ou t = 1.34). Corresponding D± 
and D 11 as a function of 4> are shown in Figs. (He) and[3Jd), 
respectively. Equation (1241) involves no free parameters, 
but reasonably agrees with the numerical results. 



IV. DISCUSSION 

We have studied diffusive light transport in a two- 
dimensional packing of monodisperse disks. We em- 
ployed ray optics to follow a light beam or photon as 
it is reflected by the disks. We used Fresnel's intensity 
reflectance with its rich dependence on the incidence an- 
gle and polarization state of the light. We note that a 
photon which moves in a grain, hits its surface with an 
incidence angle 7 less than j c — arcsin(n OMt /n.;„). Indeed 
7 < 7 C ensures that photons are not caged in grains due 
to the total internal reflection phenomena. Using a con- 
stant intensity reflectance independent of the incidence 
angle [3l| . it is not clear how the total internal reflec- 
tion phenomena influences the transport-mean-free path 
I*. Moreover, a constant intensity reflectance model does 
not take into account that n^o("f) 7^ r _>i(7). 

Our analytical estimate of D is bigger than our nu- 
merical estimate by a factor about 1.5, see Fig. O In 
writing master equations for photon transport, we do not 
consider the inequality of step lengths either in the host 
medium or in the grains. We also neglect the fact that 
the angle of refraction is not equal to the angle of in- 
cidence. But in our numerical simulations, the photon 



step length exhibits its natural distribution. Moreover, 
we strictly obey Snell's law. 

Consider a photon moving in the host medium and 
hitting a grain with an incidence angle 7. The pho- 
ton experiences an average scattering angle L (w — 
2^)F out (7)^7 = 2 (in radians) due to the reflection. 
Taking into account the Snell's law, the average an- 
gle between the incident and refracted ray is L [7 — 
a,TCsin( Ih ^ LL sin^)]F ou t(^)d'y — 0.22, where have assumed 
iii n = 1.5 and n out = 1.0. Similarly, a photon moving in 
a grain experiences an average scattering angle fj c (it — 
2j)Fin(j)dj = 2.45 due to reflection, and an aver- 
age scattering angle JJ C [arcsin( ^™^ swy) — 7]-F , m(7)<^7 = 
0.22 due to transmission. Thus, the transmission is less 
efficient than the reflection in randomizing the direc- 
tion of photons. Apparently, respecting path random- 
ization due to transmissions, our theoretical D decreases 
towards numerical one. We have also studied the distri- 
bution function of step length G(L out ). After reaching 
its pronounced maximum, G(L out ) decays exponentially, 
see Fig. 5 of Ref. [31(. Focusing on a radically different 
model, an uncorrelated random walk in a dilute pack- 
ing of point scatterers, Heidcrich et al. [41| found that a 
broader distribution of step lengths leads to a greater D. 
The interplay between the two above mentioned opposing 
impacts on D remains to be clarified. 

In two-dimensional space D = l*v m /2, where v m is 
the transport velocity of light. In a medium composed 
of spheres comparable to the light wavelength, the trans- 
port velocity differs by an order of magnitude from the 
phase velocity [42|]. The difference between the two ve- 
locities is unimportant when spheres are much larger 
than the light wavelength. Equation (|2"3"|) presents a 
"mean field" estimate of v m . Figure [3] demonstrates 
the transport-mean-free paths l*\_ = 2D±/v m and It = 
2Du/v m as a function of the packing fraction <p for the 
case nm = 1-5 and n out — 1.0. As expected, lj_ and It 
monotonically decrease as <fi increases: The photon hits 
more gains and rapidly forgets its initial direction of mo- 
tion. For packing fraction <f> « 0.64, we find a reasonable 
value I* ~ (l* ± + lf)/2 ~ 12R, cf. Refs. 0, [U, H3 • Note 
that to address photon diffusion in a three-dimensional 
system, we have taken an average over _L and || polariza- 
tions. 

To achieve a better understanding of photon diffusion 
in granular systems, there is still much to do. We aim 
at a superior model, which considers path randomiza- 
tion due to transmissions and the natural distribution of 
step lengths. Also an extension to the three-dimensional 
packing of polydisperse spheres is envisaged. Here we es- 
timate the transport-mean-free path. DWS experiments 
also reveal the dynamics of granular systems. Recent the- 
ories concerning deformations of a sphere packing, need 
an improvement as n out /ni n becomes less than 2/3 and 
the critical angle j c = arcsin(n out /rii„) deviates more 
from it/2, see Fig. 5 of Ref. [I3]. Here one must consider 
correlation between paths of a photon which moves in a 
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FIG. 3: (Color online) The transport-mean-free paths l* ± and 
Jjl (in units of the disk radius R) as a function of the packing 
fraction <j> for the case rn n = 1.5 and n ut — 1.0. Theoretical 
and Monte Carlo simulation results are denoted, respectively, 
by lines and points. 



disk and suffers many reflections [13j. Currently, we are 
studying the abrupt change of photon paths due to the 
local rearrangement of grains. Note that in another sys- 
tem, foam, bubble rearrangements cause fluctuations in 
the intensity of scattered light [26| . 

Different experiments may be performed on two- 
dimensional granular systems. One can immerse a pack- 
ing of disks in various liquids to study I* dependence on 
the refractive index n ou t- The incident electric field can 
be either perpendicular or parallel to the two-dimensional 
plane. A ray maintains its polarization state on hitting 
a disk. Quite remarkably IT > Z* , see Fig. [3] This is rea- 
sonable, since r\\ < r±: A photon sooner forgets its initial 
direction of motion as the scattering events becomes more 
probable. In another anisotropic system, nematic liquid 
crystal, the light diffuses faster along the director than 
perpendicular to the director |24| . 



Appendix A: Intensity reflectances 

Fresnel's intensity reflectance depends on the polar- 
ization state of the light. The incident electric field 
can be either perpendicular J_ or parallel || to the two- 
dimensional plane covered by the grains. For all angles 
< 7 < tt/2, 



"£*(7) = 



,(7) 



n out cos^ - m 



1 



;^»7) 2 



n out cosj + n ln 
'lout 



1 - {^^sinj) 2 - n m cos~f 



n ou tJl - (t^s««7) 2 + n ln cosj 



For all angles 7 < j c , 
rL„(7) = 



t cos7 - n out Jl - {^;sin-f) 2 



t cos7 + n out Jl - (-^sin-f) 



1 - (^^siwy) 2 + noutcos') 



(A2) 



Indeed r±> (j) = r|_> (7) = 1 if 7 > 7c 



Appendix B: The probability distribution Fi n (-y) 



A photon which moves in a grain hits its surface with 
an angle 7 < "f c , thus Fi n (pf) = if 7 > j c . To find the 
probability distribution Fi n (-j) for 7 < j c , we consider 
path of photons inside the disk, see Fig. 0Jb) . Each ray 
can be characterized by its distance s from the center 
of the disk, s = i?sin7 thus s < Rsinj c . We assume 
that the random variable s has a uniform distribution 
in the interval [0,i?sin7 c ]. The cumulative distribution 
function F c (j) = J^ F m (ip)dip is then ^(7) = Prob(s < 
i?sin7) = i? sin 7/ (i? sin 7 C ). It follows that 



Fin(l) 



dF c (rf) 

<h 



cos 7/ sin 7 C 



(Bl) 



for 7 < 7c. Further numerical simulations confirm our 
analytical results for Fi n (j) and -Fern* (7)- 



Appendix C: Functions Q°*£ and Q\ n n 



The functions Q™* and Qf™ introduced in Eq. ([TTjl. 
are the solutions of the following equations: 



Qo!n+M\m) = t^ Qt: n (a|m) + C Q^(a|m), 



Q™n+M\m) = t o ^Q o u ^a\m) + c* 



,(a|m), 



(Al) 



(CI) 



)) 



At^ o Qi%(a|m)+C ! t Q^(a|m) 

+ \e- ia {l^ Qt, n (a|m + 1) + C* 1 0^(a|m + 1 

+ ie fa (ti_>„Qfo(a|m - 1) + C'_ i g^(a|m - 1)), 

A- 1 i o ^Q^(a|m)+C, Ql^(a|m) 
+i e - ia (t ^Q^(a|m + l)+C,iQ^(a|m + l)) 

+ie ia (t ^Q^(a|m - l)+C,_ 1 Q^(a|m - 1)), 

(C2) 



set of algebraic equations whose solutions are 

( 1 - z C,o) ( 9o"n=o(a> m )+^^o < 9o™r 1 =o(" J TO ) 



r o ut (z\a,m}- 



Q l n (z\a,my 



zt o ^Qg"' =0 (a,m) + (l-^)Q™ =0 (a,m) 



A„ 



(C4) 



n ut (z\a,m)-- 



}T(z\a,m)= 



(l-zc™ )A m {a,z) + \zti->oB m {a, z) 



A 1 zt ^ i A m (a, z) + (l-zc^i)B m (a, z) 



(C5) 



\ 2 U^oQ l 2 n n {a\m) + C? Q^(a|m) 

+e- ta (\t^ Q\ n n (a\m + l)+C* 1 Q^(a|m + 1) 

+e M (AZ^ Ql';„(a|m - lHC^Q^Hm - X ) 

4™(^G£»(a|m)+c£ t Q$(a|m)) 

+i e - 2i « (i^Q™ (a]m + 2)+C t 2 Qr>l™ + 2) 

+ Je 2M (i^ Q™(a|m - 2)+C?_ 2 Q$(a|m - 2)) , 

Ovn-i(«l"») = 

+e- M (A- 1 Z ^Q™ t (a|m + l) + c™ 1 Q^(a|m+ 1)) 
+e M (A- 1 t ^Q^(«|m - l)+<_!Q5? B (a|m - 1) 



+ 2(*o-*Q5£(a|m)+c£ Qfo(a|m) 

+i e - 2 » (i^Q™ (a|m + 2)+c£ 2 Q£ n (a|ro + 2) 



+ ie ^ Q (t ^g^(a|m - 2)+C ] _ 2 Q^(a|m - 2) j , 

(C3) 



where A = Li n /L out . The above set of linear difference 
equations express Q°"* +1 and Qf l n+1 in terms of Q™f and 
Q* 7 ^ The z transform of Q n +i is simply Q(z) / ' z — Q n =a / ' z . 
Note the similarities of this rule with the Laplace trans- 
form of the time derivative of a continuous function [43l | . 
The z-transform of above difference equations leads to a 



c (z|a,m): 



l (z\a,rri)= 



(l-zc™ )D m (a, z) + X 2 zU^ E m (a, z) 



\- 2 zt ^D m {a, z) + (l-zc™%)E m {a, z) 



(C6) 



Here 



A m {a, z) = Qr„' =0 (a, m) 

+f e~ ia (i^ Q™ (a|m + 1) + C^Q° %(a\m + 1) 

+ Z ^ a (t^ Q™ n {a\m - 1) + C^Q^Hm - 1) 



B m (a,z) = Q™ {a,m) 



+ -e- M (* ^Q£ n (a|m + 1) + C^Jalm + 1 
+ ^e ia (t o ^Q° %(a\m - 1) + C^Q^Jalm - 1) 



D m (a,z) = Q° u „* =0 (a, m) 
+ze" M (At^ Ql^(a|m + 1) + c£?!Q?£(a|ro + 



+ze ra (^ Q^H TO - 1) + <\~iQ°Cn(ot\m - 1) 
t^o<#„(a|m)+c^,Q$(a|m) 

+ Z -e- 2la (t^ Q™ n (a\m + 2) + ^^(olm + 2) 
+ ^ 2M (t^oQ^Hm - 2) + c£?_ 2 QS£(a|ro - 2) 



E m (a,z) = <3™ n=0 (a,m) 
+ze- ia (\- 1 t ^ i Q° 1 %(a\m + 1)+C i g^(a|m + 1 ) 

+ze ia (\-H ^iQTn{a\m - l)+C-iQl"™H™ - I ) 



+l[t ^QZ{a\m) + c™ Q* n n (a\m] 



to-*QJ?„(a|m + 2) + C f2 QJ7 n (a|m + 2) 



and 



"4 6 



A, 



2<o 



(to-iQ^H™ - 2) + <_ 2 Qfr,(a|ro - 2) 



(1 ZC m0 )(l ZC m0 ) Z to-titi-to- V^') 



The expressions of Ql n (z\a,m) and Q°"*(,2|a:,77j) con- 
tain the sum of several terms whose inverse ^-transform 
are readily accessible: 



1 «• 



n o 



a" «• 



na <H> 



1 



1- 


5 

- z 




z 


(1 


-z) 2 ' 




1 


1- 


- az ' 




az 



{l-azf 



(C8) 



Here a is an arbitrary real number whose absolute mag- 
nitude is less than 1. 
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